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Abstract 

In a symplectic reduction framework we construct action-angle systems of canonical coordinates 
for both the hyperbolic Sutherland and the rational Ruijsenaars-Schneider-van Diejen integrable 
models associated with the C n root system. The presented dual reduction picture permits us to 
establish the action-angle duality between these many-particle systems. 
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1 Introduction 



In the theory of integrable Hamiltonian systems the construction of action-angle coordinates is of 
primary interest. The motivation comes from the fact that in these variables the equation of motion 
takes a particularly simple linearized form, whence its integration is trivial. For Liouville integrable 
systems, under certain technical conditions, the Liouville-Arnold theorem [I] guarantees the existence 
of action-angle variables, but this result is of little help in the explicit construction of these distin- 
guished canonical coordinates. The actual fact is that a bare-hand attempt for their construction 
may easily lead to non-trivial analytic subtleties. However, for certain integrable many-particle sys- 
tems defined on the real line the construction of the action-angle variables is under complete control. 
Indeed, in an ingenious paper [2] Ruijsenaars has constructed action-angle coordinates for both the 
Calogero-Moser-Sutherland (CMS) and the Ruijsenaars-Schneider (RS) models associated with the 
A n root system. One of the outcomes of his analysis is that the action variables of the non-relativistic 
models are the particle-positions of the relativistic models, and vice versa. This phenomenon goes 
under the name of action-angle duality between the CMS and the RS models. 

Though the construction of the action-angle variables has been extended to a larger class of A n - 
type particle systems (see e.g. [3], [1]), and the duality properties have been also reinterpreted in the 
context of symplectic reduction (see e.g. [5], [6], [7]), to our knowledge the explicit construction of 
action-angle variables for the non-^4 n -type CMS and RS models has not appeared in the literature 
before. The main goal of this paper is to construct action-angle variables for both the hyperbolic C n 
Sutherland and the rational C n Ruijsenaars-Schneider-van Diejen (RSvD) models. 

Let us recall that the phase space of the hyperbolic C n Sutherland model is 

V S = {(qi, . . . , q n , Ph . . . , Pn ) € M n x W 1 \ Ql > . . . > q n > 0}, (1.1) 
and the dynamics is generated by the Hamiltonian 

1 n / 2 2 \ i n 2 

H S = _ST^ 2, \ ^ / 9 9 \ 1 92 2) 

° n i<h< n Wnh 2 ^-^) sinhW^; 2^sinh 2 (2 Qc )' K '> 

where the so-called coupling parameters g and 52 are arbitrary non-zero real numbers. As is known, 
the dynamics admits a Lax representation, which in turn provides a simple solution algorithm for 
the model. Also, the non-j4 n -type Sutherland models have been successfully fitted into a convenient 
symplectic reduction framework (for details see e.g. [8], [9], [10], [11], [12], p3]). By furnishing action- 
angle variables for the hyperbolic C n Sutherland model, in this paper we complete their symplectic 
geometric understanding. 

The phase space of the rational C n RSvD model is 

V R = {(Ai, . . . , A n , 0i, . . . , 9 n ) € R n x R n I Ai > . . . > A n > 0}, (1.3) 

and the dynamics is governed by the Hamiltonian function 

n , 2\| n / a 2 \ I / /|2\| 

hs. = x> w (1 + 1) n + A?) i 1 + ot) • M 

The geometric aspects of this model is far less developed than that of the Sutherland models. Although 
the non-^4 n -type RSvD models are also known to be Liouville integrable j!4j . even the Lax represen- 
tation of their dynamics is missing. However, in our recent paper [15] we proposed a Lax matrix for 



2 



the rational C n RSvD model (|1.4p with two independent coupling parameters. In this paper we prove 
that the proposed Lax matrix (|4.2[) does have the properties one expects from Lax matrices. What 
is more important, we construct action- angle variables and also provide a solution algorithm for the 
RSvD model (|1.4p . As a by-product of our construction, we establish the aforementioned action-angle 
duality between the hyperbolic C n Sutherland and the rational C n RSvD models. 

We now briefly outline the content of the rest of the paper. Since we work out the construction of 
the action-angle coordinates in a symplectic reduction framework, section 2 is devoted to be a brief 
review on the reduction procedure. Using this machinery, in section 3 we derive the phase space of 
the hyperbolic C n Sutherland model from symplectic reduction. Though this is a standard material, 
we find it instructive to discuss the reduction picture of the Sutherland models, mainly to see the 
parallel development with the theory of the RSvD models. The presentation of the new results starts 
with section 4. First we examine the main properties of the proposed Lax matrix (|4.2|) for the C n 
RSvD model. The main new technical result of the paper is theorem 14, in which we show that the 
parametrization of the Lax matrix (|4.2[) does provide a canonical coordinate system. Having equipped 
with theorems 5 and 14, in section 5 we construct action-angle coordinates for both the C n Sutherland 
and the C n RSvD models in one stroke. By the very nature of the construction, the action-angle 
duality between these models comes for free. Also, the proposed geometric picture naturally leads to 
simple solution algorithms for both models. Incidentally, we prove the action-angle duality between 
two families of Hamiltonian systems, among which the Sutherland model (|1.2p and the RSvD model 
(|1.4p are undoubtedly the physically most interesting members. Finally, we conclude the paper with 
an appendix on some useful facts on Cauchy matrices. 

2 Preliminaries 

Take an arbitrary n E N = {1, 2, . . .} and let N = 2n. Making use of the N x N unitary matrix 



C 



On In 
In On 



e U(N), (2.1) 



we define the non-compact real reductive matrix Lie group 

G = U{n, n) = {ye GL(N, C) | y*Cy = C} (2.2) 
of dimension dim(G) = 4n 2 . Its Lie algebra has the form 

Q = u (n, n) = {Y G gl(N, C) | Y*C + CY = 0}. (2.3) 

Notice that the map 

(,): fl x fl ->R, (Y 1 ,Y 2 )^(Y 1 ,Y 2 )=tT(Y 1 Y 2 ) (2.4) 
is a well-defined Ad-invariant symmetric bilinear form on q. Since it is also non-degenerate, the map 

9 Y M- (Y, •) G * (2.5) 

is a linear isomorphism, leading to the natural identification 0* = 0. Moreover, by its Ad-invar iance, 
the bilinear form (|2.4[) naturally induces a biinvariant pseudo-Riemannian metric on the group G. 
Now let L a : y i— > ay (a G G) denote the left translations on G; then the left trivialization 

G x 0* 9 (y, ip) i y ( L i ) V G T*G (2.6) 
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permits us to identify the cotangent bundle T*G with G x g*. Moreover, due to the isomorphism 
(|2.5|) . the product manifold 

V = G x g (2.7) 

is also diffeomorphic to T*G. For convenience, in the rest of the paper we use V as an appropriate 
model of the cotangent bundle of G. Observe that the tangent spaces of V can be identified as 

T {y:Y) V = T {yX) (GxQ)^T y G®T Y Q^T y G®Q ((y,Y) eV). (2.8) 

Regarding the symplectic structure of V = T*G, the canonical 1-form G f2 1 ('P) can be written as 

ti iyX) (5y(B5Y) = (y- 1 5y,Y), (2.9) 

and for the canonical symplectic form to = — di? G fl 2 (V) we have 

oj (y! Y)(Ay © Ay, 5y © SY) = (y^Ay, 5Y) - (y^Sy, AY) + ([y' 1 Ay, y'HylY), (2.10) 

where (y, Y) G T 5 is an arbitrary point and Ay © AY, 5y © <5Y G T y G © g are arbitrary tangent vectors. 
(For more details see e.g. Proposition 4.4.2 in (16j.) 

Now, let us consider the Cartan involution O(y) = (y" 1 )* (y G G) and its fixed-point set 

K = {y G G | B(y) = y} = {y G G \ y is unitary} ^ J7(n) x J7(n). (2.11) 

Note that K is a maximal compact subgroup of G of dimension dim(i^) = 2n 2 , in which the matrix 
C (|2.ip is a central element. The Lie algebra involution 0(Y) = —Y* (Y G fl) corresponding to O 
induces the Cartan decomposition 

= «ffip (2.12) 
with the Lie subalgebra and the complementary subspace 

t = ker(# — Id) = {Y G g | Y* = — Y} and p = ker(6» + Id) = {Y G g | Y* = Y}, (2.13) 

respectively. That is, any element Y G can be uniquely decomposed as 

Y = Y + + Y_ (Y+ G t, Y_ G p). (2.14) 

Notice that dim(t) = dim(p) = 2n 2 . We mention in passing that the bilinear form fj2.4[) is negative 
definite (respectively positive definite) on t (respectively on p). 

Next, with the aid of K we introduce an isometric group action on G. Indeed, the formula 

(k L , k R ).y = k L yk-^ (y G G, (k L , k R ) G K x K) (2.15) 

defines a smooth left action of the product Lie group K x K on G. The natural lift of this action onto 
V = T*G has the form $ : (K x K) x V ->■ T 7 with 

$ (fcL)M (y, Y) = (fc £ , fo) . (y, Y) = (k L yk R \ ^Y^ 1 ) ((y, Y) G V, (k L , k R ) £ K x K). (2.16) 

For any Xl®X r G fi©C let the vector field {Xl®X r )* G X{T > ) denote the corresponding infinitesimal 
generator of <£. Obviously at each point (y, Y) G V we have 

(X L © XflJ^yj = (X L y - yX fi ) © [X R , Y] G T y G © g. (2.17) 
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Since action $ leaves the canonical 1-form (|2.9p invariant, i.e. for every group element k R ) G KxK 
we have ^j s . L = ^> the corresponding momentum map J: V — > (t © £)* takes the form 

J(y,y)(X L 0X fl ) = ^ jY) (paeXfl)^) = ((yYy-^Xz) + (-y +) Xfl>. (2.18) 
(See e.g. Theorem 4.2.10 in [IS].) Now note that the formula 

(z L © z R ,x L © x R ) m = (z L ,x L ) + (^,x fl ) (2.19) 

defines an invariant negative definite symmetric bilinear form on t © t, whence the linear isomorphism 

t(Bt3 Z L (BZ R ^ (Z L (B Z R , -) m G (t © *)* (2.20) 

allows us to make the natural identification (£©£)* = I ©I. Therefore the momentum map for action 
$ can be realized as a K x isT-equivariant map of the form 

J:7>->t©t, (y,Y)^ J(y,Y) = (yYy- 1 ) + (B(-Y + ). (2.21) 

Now, with any column vector V G satisfying the conditions V*V = N and CV + V = we 
associate the traceless N x N matrix 

£(V) = - ljv) + i(s - <7 2 )C G 1 (2.22) 

Also, let E G denote the distinguished column vector with components 

E a = -E n+a = l (aeN n = {l,...,4 (2.23) 

In order to understand the action-angle duality between the hyperbolic C n Sutherland and the rational 
C n RSvD models from a symplectic reduction picture, we wish to reduce the Hamiltonian K x iT-space 
(V,lo, <3?, J) at the very special valued 

/i = -£(E)eoetee = (eee)* (2.24) 

of the momentum map J. To perform the actual reduction, we utilize the so-called shifting trick (see 
e.g. [17j). Recall that the shifting trick relies on the co-adjoint orbit passing through — p, which in 
our case can be identified with the adjoint orbit O © {0} C £ © £ = (t © £)*, where 

O = 0{i{E)) = {£ (V) G t 1 V G C^, V*V = N, CV + V = 0}. (2.25) 

Since O = CP n_1 , for its real manifold dimension we have dim(O) = In — 2. Also, the tangent spaces 
of O can be identified as 

T p O * {ad x (p) = [X, p) | X G £} C I {p G O), (2.26) 
and the natural Kirillov-Kostant-Souriau symplectic form ui° G VL 2 (0) carried by O has the form 

u°{zAx{p)Mz(p)) = {p,[x,z]) ( P eO,x,Zet). (2.27) 

1 Note that by fixing J to the slightly more general value —£(E) © inC £ { © t depending on a new real parameter 
k £ R, we can even understand the action-angle duality between the hyperbolic BC„ Sutherland and the rational BC„ 
RSvD models with three independent coupling parameters. By working out the C„ case in this paper, we also pave the 
way for the analysis of the more involved BC n case, which we wish to publish elsewhere. 
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Following the idea of the shifting trick, we enlarge the initial phase space V by the orbit O © {0} = O, 
i.e. we consider the extended phase space 

p-t = VxQ = {( y ,Y,p)\ye G, Y G 0, p G O}, (2.28) 

and endow it with the product symplectic structure 

u cxt = uj + u° . (2.29) 

The natural diagonal action of K x K on V cxt is given by <£ cxt : (K x K) x V cxt -4 T 3 ^, where 

*(£,fc«)(f' y ^) = (^,^,).(y,y,p) = (kLyk^^knYk^^Lpkl 1 ), (2.30) 

and the corresponding if x i^-equivariant momentum map J ext : T 76 ^ —>•£©{ takes the form 

J ext (y,F,p) = ((yry- 1 ) + + p) © (-Y+). (2.31) 

Notice, however, that tr((yYy~ 1 ) + + p)+tr(— Y + ) = 0. Therefore, upon introducing the Lie subalgebra 

s(i © t) = {X L © X R G fi © 1 1 tr(X L ) + tr(X fl ) = 0} < « © t, (2.32) 

we can write J ext (y,Y,p) G s(tffi£) (V(y,Y,p) G T 70 ^); that is, J ext is actually an s(t © C)-valued map. 

The shifting trick, which basically states that the reduced phase space V // ^{K x K) is symplecto- 
morphic to the symplectic quotient "P ext //q(K x K), is only a convenient technical tool in our present 
work. In general this approach does not make the reduction simpler, only relegates the differential- 
topological difficulties into the orbit part of V . However, our experience with the Calogero-Moser- 
Sutherland-type many-particle systems convinces us that the shifting trick does provide a cleaner and 
shorter derivation of these integrable systems in the symplectic reduction framework (see e.g [12J, 

Finally, recall that the Hamiltonian vector field X h G X(V ext ) associated with an arbitrary smooth 
function H G C°°(V ext ) is (uniquely) determined by the condition 

dH = X H ju ext , (2.33) 

and the Poisson bracket of any pair of functions F, H G C°°(V ext ) is defined as 

{F, H} cxt = u; cxt (X F , X H ) = X H [F}. (2.34) 

3 The phase space of the hyperbolic C n Sutherland model 

The first step of the Marsden-Weinstein reduction is to solve the momentum map constraint 

J cxt (y,Y,p) = (3.1) 

for (y, Y, p) G V cxt = G x g x O; that is, we have to understand the properties of the level set 

Z = (J cxt )-HM) CV cxt . (3.2) 

Our first goal in this section is to solve the constraint (|3.1|) by diagonalizing the group elements y G G. 
However, to make the notion of diagonalization precise, we need some more Lie theoretic facts, which 
we summarize in the following subsection. For a general reference on the related group theoretic issues 
we recommend 1181. 
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3.1 Group theoretic background 

For any real n-tuple q = (q%, . . . , q n ) G IR n we define q = diag(gi, . . . , q n ) and Q = diag(q, — q) G p. 
Clearly the set of diagonal matrices 

a = {Q ep\ q eR n } (3.3) 

forms a maximal Abelian subspace in p. If denotes the set of the off-diagonal elements of p, then 
we have the orthogonal decomposition p = o © cr 1- . Note that the centralizer of a inside K is the 
Abelian Lie group 

M = Z K (a) = {diag(e ix , e ix ) G K \ X G R n } < K (3.4) 

with Lie algebra m = {diag(i%, i^) G t \ x G M. n } < i. If m -1 denotes the set of the off-diagonal elements 
of 6, we can write 6 = m © m^. Since the subspace © a 1 - consists of the off-diagonal elements of g, 
the restricted operator 

ad Q = ad Q | m ± ffia ± G Ql(m x © cr 1 ) (3.5) 
is well-defined for any q G W 1 , and its spectrum has the form 

cr(adg) = {q a - qb, ±(q a + Ob), ± 2 Qc \ a, b, c G N n , a / b}. (3.6) 

The regular part of a is defined as a reg = {Q G o | adg is invertible}, and the standard Weyl chamber 

c = {Q G a rcg | qx > . . . > q n > 0} (3.7) 

is an appropriate connected component of a reg . Without any further notice, in the rest of the paper 
we frequently identify c with the subset {q = (q\, . . . , q n ) G M. n \ q\ > . . . > q n > 0} C M n . 

At the Lie algebra level it is a crucial fact that the elements of p (|2. 13[) can be conjugated into a 
by appropriate elements of K. More precisely, the map 

o x K 3 (Q, k) ^ kQk- 1 G p (3.8) 

is well-defined and onto. The regular part of p defined by p reg = {kQk^ 1 | Q G c, k G K} is dense and 
open in p. Moreover, the map 

c x K/M 3 (Q, kM) ^ kQkT 1 G p rcg (3.9) 

is a diffeomorpism, allowing us to make the identification p reg = tx K/M. 

To proceed further we need some basic facts from the canonical forms of the Lie group elements, 
too. First, recall that the global Cartan decomposition (polar decomposition) 

exp(p) xK3 (e Y , k) ^e Y k£G (3.10) 

is a diffeomorphism. Combining this fact with (|3.8p . it is obvious that any group element y G G can 
be written as 

y = k L e Q k^ (3.11) 

with some kL,kR G K and Q G a. Upon setting A = exp(o), the factorization ([3. lip of the elements 
of G is usually called the KAK decomposition. 

Let us now define the regular part of the Abelian Lie group A by = {e^ G A \ Q G a reg }, and 
also let G reg = KA reg K. It is known that G reg is a dense and open submanifold of G. Moreover, if 
M* denotes the diagonal embedding of M (|3.4p into K x K, i.e M* = {(m,m) G K x K \ m G M}, 
then the map 

cx(Kx K)/M* 3 (Q, (k L , k R )M*) ^ k L e Q k^ G G reg (3.12) 
is a diffeomorphism. Thus the identification G Teg = c x (K x K)/M* is immediate. 
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3.2 Parametrization of £ induced by the KAK decomposition 

Having equipped with the necessary group theoretic facts, now we are in a position to derive both 
the hyperbolic C n Sutherland and the rational C n RSvD models from a uniform symplectic reduction 
framework. To derive the Sutherland system (|1.2p defined on (jl.ip . one of the key ingredients is the 
Lax operator 

L:V s ^p, (q, P )^L(q,p)=P-smh( a d Q )- 1 Z(E), (3.13) 

where we employ the notation P = diag(p, — p) S o with p = diag(pi, . . . ,p n ), as we stipulated in the 
previous subsection. Notice that the matrix entries of L have the form 

L a ,b = —L n -\- a ,n+b = . , / s" j ^a,n+b = ~L n +a,b = . , / ; C > (3-14) 

smh(g a - q b ) smh(g a + q h ) 

Lc,n+c — L n J rCC — — j— r, L cc — L'n+c.n+c — Pci (3.15) 

smn(2g c J 

where a, b, c € N n and a ^ b. Now, making use of the Lax operator L, the KAK decomposition of G 
leads to the following characterization of the points of the level set £o (|3.2p . 

Proposition 1. For each (y,Y,p) € £o there exist some q G c, p E M. n and til,tir € K such that 

y = r lL e Q ri R 1 , Y = r) R L(q,p)rf^ , p = n^E)^ 1 . (3.16) 



Proof. Take an arbitrary (y,Y,p) £ £o; then by the KAK decomposition (|3.11|) we have the factor- 
ization 

V = ^fj' 1 (3.17) 

with some f]L^R € K and Q = diag(q, — q) € o satisfying q% > . . . > q n > 0. Also, due to ()2.25j) . 
we can write p = £(V) with some V G satisfying V*V = N and CV + V = 0. Plugging this 
parametrization into J ext (|2.3ip . we conclude that Y + = 0, i.e. Y = Y_ € p, together with 

= (yyy- 1 ) + + £(V) = fiL (sinh(ad Q )(^ 1 y_^,) + ^V)) rfc 1 . (3.18) 

Upon introducing Y = fj^Y-fjR G p and V = f\J^V G C^, the above equation can be rewritten as 

sinh(adQ)y = -£(F). (3.19) 

Componentwise, for any k,l € Ntv we have 

sinh(g fc - qi)Y k! i = ig(S k ,i - V k Vi) + i(g 2 - g)C kjh (3.20) 

where it is understood that q n + a = —q a (Va G N n ). 

Take an arbitrary a £ N n ; then with k = I = awe get = i<7(l — | | 2 ) , therefore = —V n+a = e 1Xa 
with some Xa € M. With k = a and / = n + a we obtain sinh(2(7 a )y ajra+a = ig 2 ^ 0, whence q a ^ also 
follows. 

Next, let a, 6 £ N n , a ^ b; then with k = a and I = b we have sinh((7 a — (ft)y aj & = — igV a Vb ^ 0, 
which implies g a 7^ qv At this point we see that actually 56 c. Moreover, upon introducing the group 
element m = diag(e lx , e lx ) € M, the application of the linear operator Ad m -i on (|3.19p yields 

sinh(ad Q )Ad m -iy = -£(E) € m ± , (3.21) 
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therefore Ad m -iY = P — sinh(adQ) _1 ^(-B) = L(q,p) with some p E M. n . Finally, by letting t]l = fj^m 
and r]R = fjRm, the parametrization (|3.16p immediately follows. ■ 
The characterization of the points of the level set £o given in proposition 1 permits us to introduce 
a parametrization naturally induced by the KAK decomposition of the group elements. For, we need 
the diagonal embedding of U{\) into K x K; i.e. the Lie subgroup 

?7(1)* = {(e ix l N , e lx l N ) £K xK\ X £^} = U(l), (3.22) 

and the corresponding Lie algebra 

u(l)* = {(i X lN, ixliv) G t 1 1 x e K} = u(l). (3.23) 

Also, let us define the smooth product manifold 

M s = V s x (if x K)/U{1)*. (3.24) 

Since 17(1)* is a normal subgroup of If x K, the coset space (K x K)/U(l)* is a (real) Lie group in a 
natural manner of dimension 4n 2 — 1. Note also that dim.(A4 s ) = 4n 2 + 2n — 1. 

Lemma 2. T/ie map 

T 5 : A4 s ^P ext , (?,p,(77L,%)^(l).) ^ fa^SaS vrL(q,p)Vr\ VlZWvZ 1 ) (3-25) 

is an injective immersion with image T S (M S ) = So- 
Proof. It is obvious that T is a well-defined smooth map. Our next goal is to determine the image 
of M s under the map I s . To this end let q E c, p E W 1 and t]l,7]r £ if be arbitrary elements and 
consider the matrices 

FvSs'eG, Y = VrL(q,p)Vr 1 6P, P = %£(£)% 1 £ O. (3.26) 

Notice that the relationship 

yYy- 1 = ry L (cosh(adQ)L(g,p) + sinh(adQ)L(g,p)) r/J 1 (3.27) 

together with ()3.13[) immediately yields 

(yyy- 1 )+ + p = ? ?L (sinh(adQ)L(a,p))^ 1 + ^{E)^ 1 = 0. (3.28) 

Since the equation Y + = also holds, from the definition of the momentum map J ext (|2.3ip we see at 
once that J ext (y,Y, p) = 0. In other words, we have the inclusion relation T S (A4 S ) C £q. Recalling 
proposition 1 we conclude that T S (A4 S ) = £o- 
To show that T s is injective, suppose that 

T 5 (q,P,(VL,r]R)U(l)*) = T S (q,p,(rj L ,fj R )U(l)*) (3.29) 

with some q,q £ c, p, p £ M n and rjL,fjL,r]R,fjR £ K. The G-component of the above equation reads 
as 

(VLe Q rjl 1 )(rjLVR 1 ) = VLe Q v R 1 = VL^fjl 1 = feeS" 1 )^^ 1 )- (3.30) 
Remembering the uniqueness of the global Cartan decomposition (|3.10p . the relations 

i]Le Q i] r [ 1 = fibril 1 and = Vl^r 1 (3.31) 
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are immediate. Moreover, from the regularity assumption q, q G c and from the identification (|3.9p we 
see that q = q and t)l = rum with some m G M. Notice that the relationship fj R = r] R m also follows. 
The 0-component of (|3.29p has the form 

VL^E)r,l 1 = ij L ^E)ri£ l , (3.32) 

hence £(E) = m£(E)m _1 = £(mE). By ([2321) we can write EE* = (mE)(mE)* , which immediately 
implies that m = e lx ljv with some x € M. Therefore (r]L, 7/#)t7(l)* = (f]L,f]R)U(l)*. 
Finally, the g-component of (|3.29[) spells out as 

r lRL(q,p)VR 1 = '>1rL{QjP)Vr 1 ■ (3-33) 

Utilizing the above observations, this equation leads to L(q,p) = L(q,p). By simply taking its diagonal 
part, we obtain P = P, i.e. p = p, completing the verification of the injectivity of T s . 

Next, we show that T s is an immersion, i.e. its derivative is injective at each point of M. s . For, 
take an arbitrary point x = (q,p, (r]L,r]R)U(l)^) G Ai s . Due to the natural identifications 

T q c^R n , T p R n = M. n , T e ((KxK)/U(l)*)^(t(Bt)/u(l)*, (3.34) 

a generic tangent vector v at the point x has the form 

v = 5q (B 5p (B (vl,Vr)(Xl ® X R )u(l% G T X M S (3.35) 

with some 5q,5p £ M. n and Xl,X r G t. Notice also that the action of the derivative 

(dT s ) x : T X M S T rS{x) V ext (3.36) 

on v is given by the formula 

(dX 5 )^ = vUe^SQ + X L e® - Xr)^ 1 ® m {8L + [X R , L{q,p)])rf^ Vl [X l , {(E)}^ 1 , (3.37) 

where 5L = 5P + [SQ, cosh(adg) sinh(adQ) _2 £(£)]. 

Now, to find the kernel of the linear map (dT s ) x , suppose that v G ker((dT s ) x ), i.e. (dT s ) x v = 0. 
Taking the G-component of (dT s ) x v, for any k,l G Nat, k 7^ /, we obtain 

(X L )k,ie gi =e q *(X R ) k j, (3.38) 

with the convention <7 n +a = — q a (Va G N n ). Combining (^38]) with X£ + X L = and X* R + Z K = 0, 
it follows that 

Bmh(q k - qi )(X R ) ktl = 0. (3.39) 

Since q is regular, we see that both Xl and X R are diagonal matrices, i.e. Xl,X r G m. Plugging this 
observation back into the G-component of (dT s ) x v, we obtain 5Q + Xl — X R = 0. Inspecting its real 
part we get 5Q = 0, i.e. 5q = 0, meanwhile the imaginary part yields Xl = X R . 

Remembering the defining formula of £(E) (|2.22p . from the vanishing of the O-component of the 
derivative (|3.37p we see that 

= [X L ,S(E)\ = ig[X L , EE*} = ig {(X L E)E* + E(X L E)*) . (3.40) 

Therefore, for any a, b G N n we can write 

= {X L E) a E b + E a JX~L~E) b = (X L )a,a + {XL)b,b- ( 3 - 41 ) 

It follows that Xl = X R = ix^-N with some x € M, i.e. Xl © X R G u(l)*. Finally, by inspecting the 
g-component of (|3.37p . we obtain 8P = 0, i.e. 5p = 0. Thus v = 0, meaning that the kernel of the 
derivative operator is trivial. ■ 
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3.3 Identification of the reduced phase space 

In this subsection we wish to identify the phase space of the hyperbolic C n Sutherland model in the 
proposed symplectic reduction picture. As a first step, we examine the differential geometric properties 
of the closed level set £o (|3.2p . 

Lemma 3. The zero element o/s(t©fi) is a regular value of the momentum map J ext : J> cxt — y s(£©{). 
Proof. Pick an arbitrary x = (y,Y,p) G •p ext . Working out the derivative operator 

(dJ cxt ) x : T x V cxt -> Tj^ {x) s(i © t) ^ s{t © t) (3.42) 

of the momentum map J ext (|2.3ip at point x, for each tangent vector 5y © SY © Sp 6 T x V ext we find 
the formula 

(dj cxt ) x (5y © © M = (G/(«yy - [r, y-^yDy -1 )* + s P ) © (-(«-)+) . (3.43) 

Now, we have to verify that the derivative (|3.42|) is onto at each point x G £o- Since J ext is equivariant, 
by lemma 2 it suffices to show that (|3.42p is onto at each x = (e^, L(q,p), £(E)), where q G c and 
p G MJ 1 . For this purpose, choose an arbitrary element W © W G s(t © t) and decompose it as 

W ®W = {A + ml N )®(B -\nl N ), (3.44) 

where are traceless matrices and k G R. Upon defining the anti-Hermitian tridiagonal matrix 

T € M nxn with the only non-zero entries 

a A + B 

T aA+1 = -T a+lA = C ' C 9 c ' c £1 (a G N n _i), (3.45) 

c=l ^ 5 

we see that the N x N matrix X = diag(T, T) belongs to the subalgebra Therefore X generates a 
well-defined tangent vector 5p = [X, £(E)] G T^rmO. Moreover, let 5y = and consider the tangent 
vector 5Y = (5Y)+ + (8Y)~ G determined uniquely by the requirements 

(5Y)+ = -B + iKl N = -W and ((<Jr)_) dia€ = 0, (3.46) 

together with 

((5F)_) ofr _diag = sinh(adQ)" 1 (^ off _dia g - (<fy>)off-diag - cosh(adQ)((<5y) + )off-diag) G a -1 - (3.47) 

By applying the derivative operator on the distinguished tangent vector 5y © SY © 5p defined by the 
above conditions, a straightforward calculation based on (|3.43|) yields 

(dJ cxt ) x (5y ®5Y ®Sp) = W®W. (3.48) 

In other words, the derivative is onto. ■ 
Corollary 4. There is a unique smooth manifold structure on £o 13. 1|) such that the pair (£o,to) 
with the natural inclusion lq : £o 'P ext is an embedded submanifold of , p ext . Moreover, the pair 
(A4 S ,T S ) provides an equivalent model for (£q,lq), i.e. £o — M S ■ 

Proof. As is known (see e.g. Proposition 1.6.18 in |16j). the previous lemma guarantees that the level 
set £o C "P ext is an embedded submanifold of dimension 

dim(£ ) = dim(7> ext ) - dim(s(« © I)) = An 2 + 2n - 1 = dim(A4 5 ). (3.49) 
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From lemma 2 we also see that T is a smooth bijective immersion from A4 S onto £q, therefore they 
are necessarily diffeomorphic. ■ 
Utilizing the model (A4 S , T s ) of the level set £q, in the remaining part of this section we complete 
the reduction of 'P ext at the zero value of the momentum map J ext . Note that the stabilizer subgroup 
of E s(i ffi t) is the full group K x K, and the (residual) K x K action on A4 S takes the form 



(k L ,k R ) . (q,p, (r)L,r) R )U(l)*) = (q,P, (kLVL,k R r] R )U(l)*) , 



(3.50) 



for any q G c, p G M. n and kL,k R ,i]L,rj R G K. Thus, it is obvious that the orbit space A4 S /(K x K) 
can be identified with the base manifold of the trivial fiber bundle 



(3.51) 



vr S : M S ^V S , (q,p, (r, L , Vr)U(1)*) -> (q,p), 

therefore the reduced symplectic manifold can be identified as 

V cxt // (K x K) = M S /(K x K) = V s . (3.52) 

Recall also that the reduced symplectic structure ui s G Q 2 (V S ) is uniquely determined by the condition 

(n 3 )*^ = (T 5 )*6j ext . (3.53) 

Working out the above pull-backs, the reduced symplectic form can be determined in a straightforward 
manner. 

Theorem 5. Up to a multiplicative constant, the globally defined coordinate functions q a , p a (a G N n ) 
provide a Darboux system on the reduced manifold V s . More precisely, the reduced symplectic form 
can be written as u s = 2 Yla=i ^Qa A dp a . 

Proof. Take an arbitrary r = (q,p) G V s and let x = (q,p, (In, ljv)^(l)*) £ -M s . Observe that the 
point x projects onto r, i.e. ir s (x) = r. Also, we have 



dq c 



d_ 

dq c 



and (d7r 

dp c 



_d_ 

dp c 



(3.54) 



for every c G N n . If G gl(N,C) (k,l G Nat) denotes the elementary matrix with entries {e.k,l)k',V 
Skk'&ll'i then by ([3.37|) we can write 



dq c 



d 



e^W c © [W c , Z]eO and {dT b ) x — 

dp c 



© W c © 0, 



(3.55) 



where 



W c = e CiC — e n+Cin+c G o and Z = cosh(adg) sinh(a~dg) 2 £,(E) G t. (3.56) 
Now, remembering the defining formulae of w ext (|2.29p . from (|3.53|) it is immediate that 



LU,. 



(- 


d 


\dq a 





u$ (x) (e Q W a © [W a , Z] © 0, © W b © 0) = (W a , W b ) = 2S a>b , (3.57) 



for any a, b G N n . The remaining components of u s come along the same line. ■ 
The above theorem is well-known in the literature (see e.g. |11| . |12|). After developing an analogous 
theorem for the rational C n RSvD model in section 4, its importance will be transparent in section 5. 
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4 The phase space of the rational C n RSvD model 

Starting with this section we present our new results on the C n -type rational RSvD model with two 
independent coupling parameters. Our first goal is to derive the natural phase space of the RSvD 
model from symplectic reduction. As the initial step of the reduction, we have to explore the properties 
of the level set £o f|3.2|) . Contrary to the Sutherland case, now we are solving the momentum map 
constraint (|3.ip by diagonalizing the Lie algebra component of the extended phase space (|2.28p . 
First, for each a G N n let us consider the complex- valued rational function 

(<¥a) 

Remembering the manifold V R introduced in (|1.3|) . we also define the matrix- valued function 

A: V R ->gl(N,C), (\,9)^A(\,9), (4.2) 
where the matrix entries lying in the diagonal n x n blocks are given by the formulae 

A a , b {\e) = e e ^\z a (\)z h {\)\\ 219 , (4.3) 

215 + X a - x b 

A n+a n+b (X, 0) = e -'-'» 5Eg(A) 2i, 

\z a (\)z b (\)\2 2,g-\ a + \ b 

meanwhile the matrix entries belonging to the off-diagonal n x n blocks have the form 

A a , n+b (X,9) = A n+b , a (\,9) = e d «~ e >z b (\)\z a (\)z b (\r 1 \ l 2 . 2l f + l{9 ~ 92) 5 atb , (4.5) 

2ig + X a + X b ig + X a 

for any a, b G N n . In our paper [15] we introduced (|4.2|) as a natural candidate for the Lax matrix of 
the rational C n RSvD model, and we are now ready to prove this claim. As one of the key objects of 
our discussion, in the following subsection we examine the main properties of A in detail. The Cauchy 
matrices play a prominent role in our analysis, so the reader may find it useful to consult the appendix 
for the necessary background material. 

4.1 The properties of the Lax matrix A 

To facilitate the comparison with our earlier work on the structure of A, we borrow some notations 

from [TS]. First, for each A G 1" let x 6 denote the column vector with the purely imaginary 
components 

x a = -x n+a = ^ € iR (a € N n ). (4.6) 
Upon introducing the real parameter 

e = 1 - G K, (4.7) 

the function z a (|4.ip can be easily related to the rational function w a (|A.10p naturally appearing in 
the study of the standard Cauchy matrices (|A.2|) . In fact z a is a deformation of w a , as can be seen 
from the relationship 



*«(A) = -w a (x) ( 1 - i + 2x ) (a G N n ). (4.8) 
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Next, for any (A, 9) G V R we define the column vector J-(X, 9) G with components 



T a {\9) = e^\z a {\)\2 and T n+a {\9) = z a {X)\z a (X)\^ , (4.9) 
where a G N n . Then the entries of A (|4.2p can be succinctly written as 

A Kl = T f' + £ ° k ' 1 (k,lEN N ). (4.10) 
1 + x k - Xi 

Giving a glance at (|4.10p . it is obvious that A(X,9) is a Hermitian matrix. Our next goal is to show 
that A(X,9) belongs to the Lie group U(n,n) (|2.2p . 

Proposition 6. For every (A, 9) G V R we have A(X,9)CA(X,9) = C, i.e. A(X,9) G U(n,n). 
Proof. Pick an arbitrary (A, 9) G V R . Let us introduce the column vector u(x) G with components 

M*) = r^— ( k G N ^)' ( 4 - n ) 

1 + 2x k 

and the JVxJV diagonal matrices 

D = diag(J r i, . . . , Tn) and V{x) = diag(iii(x), . . . , un(x)). (4-12) 

Note that Cu{x) = u(—x) and CT>(x)C = T>(—x). Recalling (|4,10p we can write 

A = DC(x)D* + eV(x)C, (4.13) 

where C(x) stands for the standard N x N Cauchy matrix (IA.2P associated with the vector x G (iM)^ 
(|4.6|) . Thus it is clear that 

AC A -C = D(C(x)D*CDC(x)D*CD + eC(x)CV(x)CD*CD 
+ eV(x)C{x)D*CD + e 2 V(x) 2 - l^D^C. 

However, remembering (|4.9p . (|4.8p and ()A.6p . we can write 

D*CD = Cdiag(zi, ...,z n ,z 1 ,...,z n ) = CW{x){eV{x) - 1 N ). (4.15) 

Therefore, plugging this relationship into (|4.14p and applying (|A.12p . we obtain 

ACA-C = eD{l N + C{x)CW(x))V(x)(l N + C{x)CW{x)){eV(x) - l^D^C. (4.16) 

Now let us notice that 

C{x)CW{x)V{x) + V{x)C{x)CW{x) = 2u{x)u{x)*W{-x)C, (4.17) 

from where it is immediate that 

(1 N + C{x)CW{x))V{x){l N + C{x)CW(x)) = 2C (u{-x) + C(-x)W(x)u{x)) u(x)*W{-x)C. (4.18) 

However, from (|A.llj) we see that for each k G Nat we can write 

N 

(u(-x) +C(-x)W(x)u(x)) k = u k (-x) + ^2C k j(-x)wj(x)uj(x) = 0, (4.19) 
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which immediately leads to the conclusion AC A — C = 0, i.e. AC A = C . ■ 
To proceed further, let us note that the image of p (|2.13|) under the exponential map can be 
identified with the positive definite elements of U(n,n), i.e. we have 

exp(p) = {y e U(n, n)\y> 0}. (4.20) 

Keeping in mind this fact, the proof of the following lemma is quite straightforward. 

Lemma 7. For every (A, 0) G V R the Hermitian matrix A(X,6) belongs to exp(p). 

Proof. Take an arbitrary point (A, 9) G V R and keep it fixed. Note that when e = 0, the lemma is 
trivial. Indeed, in this special case A (|4.10p is a Cauchy-type matrix, thus a direct application of the 
determinant formula (|A.3j) immediately shows that all its leading principal minors are positive. 

In order to verify the lemma for arbitrary values of the coupling parameters g and g2 , let us notice 
that the dependence of A (|4.10|) on the real parameter e (j4.7[) is continuous. Moreover, according 
to proposition 6, if we change this parameter continuously from to an arbitrary given value, then 
during the course of the deformation the matrix A remains invertible and Hermitian. Therefore, by 
continuity, the eigenvalues cannot pass through zero, i.e., they remain positive. It is now evident that 
for all e the Hermitian matrix A is a positive definite element of U(n,n), whence by the identification 
(|4.20p the proof is complete. ■ 

An immediate consequence of the above lemma is that A(X, 9) has a (unique) positive square root 

K(\,0) = A(\,0)%, (4.21) 

which also belongs to exp(p). In particular, the column vector 

V(A, 0) = A(X, 0)~^J*(A, 6) G C N (4.22) 

is well-defined for all (A, 9) G V R . The importance of the following algebraic properties of V will be 
clear in the next subsection. 

Proposition 8. For any (A, 9) G V R we have V(A, 9)* V(A, 9) = N and CV(A, 9) + V(A, 9) = 0. 
Proof. Recalling (|4.8p , (|4.9p and (|4.10p , the definition of V (|4.22p leads to the following expansion 



N 



V*V = J 7 * A 



w k (-x)wi(x) 



Xk + Xl 



k,l=l 

N i \ / \ / i \ N 



£ 



+ £ 2 



+ \l-2x k l + 2xj ^ i l + 2x l ] (! - , ' ) 
y> w k (-x)wi(x) 1^ wi{x) 



Due to (|A.5p and ()A.7p . for the first sum we can write 

N , s , s N N , s N 

1 — Xh + Xl 1 + Xl — Xh 

k,l=l K 1 k=l 1=1 1 K k=l 
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Making use of (|A,11|) . the coefficient of e can be rewritten as 

A w k (-x)wi(x) ( 1 1 \ A Wj(x) _ 2 ^ w k (-x) ^ wi(x) _ 

^ l-x k +xi I 1 - 2x k I + 2xi J 1 + 2x\ ^ 1 - 2x k ^ 1 + 2xi ' K ' 

k,l = l K < l \ « 1/ l=1 I k=1 1=1 

Finally, using (|A.5p and (jA.llj) . for the coefficient of e 2 we obtain 



E w k {-x)wi{x) 1 wi(x) 
, i=i l-x k + Xl (1 - 2x k ){l + 2 Xl ) ^ (1 + 2x ; ) 2 " 

= / yv w k (-x) 1 

~ ^ 1 + 2x, \ (1 



(4.26) 



^l + 2 Xl \j^ i {l-x k + xi){l-2x k ) l + 2 Xl 



Putting the above formulae together we end up with V*V = N. 

Next, notice that CV + V = A~ a (ACT + F) . However, recalling identities (|A.5|) and (|A.11|) . for 

any k E Njv we can write 



+£ y 

, i ' x j + x k (1 - Xj + - 2xj) l + 2x k 

therefore the relationship CV + V = also follows. 



4.2 Parametrization of the level set £ 

At this point we are in a position to introduce an appropriate parametrization of the level set £o 
(|3.2p naturally induced by the diagonalization of the Lie algebra part of "P ext (|2.28p . Indeed, take an 
arbitrary A = (Ai, . . . , A n ) E W 1 and introduce the N x N diagonal matrix 

£(A) = diag(A, -A) = diag(A a , . . . , A„, -Ai, . . . , -A n ) E a, (4.28) 

then the points of £o can be characterized as follows. 

Proposition 9. For each (y,Y,p) E -Co there exist some (A, 0) E V R and r]L,7]ji E K such that 

y = r, L A(X,0)hR, Y = r ]R C(\)i lR \ p = r ]L ((V(X,8))rjl 1 . (4.29) 



Proof. Take an arbitrary (y,Y,p) E £o- Since J ext (y,Y,p) = 0, from (j2.31[) we see that Y + = 0, i.e. 
Y = y_ E p. Therefore, by (|3.8p . we can write 

Y = m^WVn 1 (4.30) 

with some tjr E K and A = (Ai, . . . , A n ) E M n satisfying Ai > . . . > A n > 0. 

Also, due to the global Cartan decomposition ()3.10p . the group element y can be uniquely factorized 
as y = y~y+, where y_ E exp(p) and y + E K. Upon introducing tjl = y+rjR E K, let us observe that 
fh^Vm = 1 lL 1 y- r lL e exp(p), therefore 

j, = 7 ?L e A ^ 1 (4.31) 
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with some A £ p. 

Plugging (|4.3U|) and (|4.31|) into J cxt (y,Y,p) = 0, and recalling (|2.25|) . we also get 

= (yYy) + + p = VL (sinh(ad A )£(A)) r ?L 1 + £(V), (4.32) 

where p = £(V) with some complex column vector V G satisfying V*y = AT and CV + V = 0. 
Note that equation (|4,32p is equivalent to 

sinh(ad A )/:(A) = -^V), (4.33) 

which immediately translates into 

2kje 2A + £(A)e 2A - e 2A £(A) = 2i 5 (e A 7 ? ~ 1 y)(e A ^V)* + 2i( 5 - g 2 )C. (4.34) 

Fortunately the last two equations coincide with the pivotal relations that permitted us to understand 
the scattering theory of the hyperbolic C n Sutherland model. Indeed, compare the above relations 
with equations (28) and (31) in |15j . Notice also that these relations are the natural analogs of the 
commutation relations in Ruijsenaars' treatment on the action-angle duality between the hyperbolic 
Sutherland and the rational Ruijsenaars-Schneider models of type A n (see equation (2.4) in [2]). 
Without repeating the arguments presented in section 3 of [15] , we simply quote the outcome of the 
analysis of equation ()4.34p . First, due to lemma 1 of [15], £(A) is a regular element of p, therefore 
Ai > . . . > A n > 0, i.e. A G c. Second, by lemma 2 of [15], we can write e 2A = A(\,0) with some 
9 G M. n , whence by equation (|4.3ip we obtain 

y = VL A(X,e^r ]R \ (4.35) 

Third, the analysis yields the relation e A n L " 1 V = T(X,6), too, therefore V = r}LA(X,6)~2j r (X,6), i.e. 

P = aVLV(X,e)) =VLtO>{X,6))vI 1 - (4-36) 

Having a look at on equations (|4,30p . (|4.35p and (|4.36p . the proposition follows. ■ 
Obviously the above proposition is the complete analog of proposition 1. Now let us introduce the 
manifold 

M R = V R x (K x K)/U(l)*. (4.37) 

Developing the theory parallel to the previous section, the following lemma can be seen as the natural 
analog of lemma 2. 
Lemma 10. The map 

T R . M R ^ £()5 (A; 0} (ri L , m )U(l)J ^ ( V lA(X, e^n^^RCiX^R^vUiViX, e))^ 1 ) (4.38) 

is a diffeomorphism, i.e. the pair (A4 R , T R ) provides an appropriate equivalent model for the embedded 
submanifold (£o^o)- 

Proof. First, we check that is a well-defined map. For, take an arbitrary (A, 0, (t]l, 77k)£7(1)*) G 
M R . As we have discussed after lemma 7, the square root of A(X,6) does belong to exp(p), whence 
we have 

y = i lL A(X,6) l 2r ]R 1 eG. (4.39) 

2 To be precise, in [13] this fact is proved under the mild assumption gi ^ 2g. Because of the quadratic appearance 
of the coupling parameters in the Hamiltonian systems of our interest (see equations (|1.2[) and (TTJ4])), this technical 
condition does not mean any restriction on our results. 
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Because of proposition 8, the column vector V(A, 9) ()4.22p generates a well-defined element of the orbit 
O flZggD , therefore p = r?L£(V(A, fl))^ 1 G O. Also, note that Y = 7] R C{\)r]^ G p. To proceed further 
we evaluate J ext (|2.31|) at the point (y, 1", p) G "P cxt with the above components. Notice that 

(yYy- 1 ) + + p= ±ri L (A(\,e)*£(\)A{\,e)-* - A(\,e)-*C{\)A{\,0)* + 2£(V(A,0)))f£ 1 . (4.40) 
However, due to ([232]) . (|4T22j) and (jiTTOj) . we can write 

A%CA~% _ A~*£J& + 2£(V) = -2K7./H + [(215)- 1 /:, .4] - J\T - eC) >H = 0, (4.41) 

therefore {yYy~ 1 ) + + p = immediately follows. Combining this fact with the obvious relation 1+ = 0, 
we see that J ext (y,Y, p) = 0, i.e. is well-defined. 

Obviously is a smooth map. Moreover, as we can infer from proposition 9, it is surjective. To 
complete the proof of the lemma, let us notice that T R acts between manifolds of the same dimensions. 
Thus it is enough to show that T R is an injective immersion. 

Let (A, 9, (r]L,r]ji)U(l)^) and (X,9, (f]L, fjji)U(l)*) be some points of M. R and suppose that 

T R (\,0,(VL,VR)UQ)*) = T R (X,9,(fj L ,fj R )U(l),). (4.42) 

The g-component of the above equation has the form 

ri R C(X) V ^=fj R C(X)f}^. (4.43) 

Since A, A G c, from (|3.9|) we obtain £(A) = £(A), i.e. A = A. Moreover, we can write fj R = r] R m with 
some unique group element m = diag(mi, . . . , mjv) G M (|3.4p . 
Next, notice that the G-component of (|4.42p yields 

(^(A,^^ 1 )^ 1 ) = r, L A(X,9)-2r,^ = fj L A(X,9)k^ = (fjLACxJ^fj^XfjL?)^). (4.44) 

By the uniqueness of the global Cartan decomposition (|3.10p . we obtain fji, = rjim together with the 
relationship 

A(X, 0)3 = mA(X, e)im _1 . (4.45) 
Now, by (|4.10p . the matrix equation A(X, 9) = mA(X,9)m~ 1 immediately leads to the relations 

T k (X,9)J^X^ = m k T k (Xj)J^xJ)mY 1 {k,leN N ). (4.46) 

Recalling the definition (|4,9p . with k = I = a (a G N n ) we get e 26a = e 29a , whence 9 = 9. Plugging 
this back into (|4.46p . for all k, I G Nat we obtain m k = mi, therefore m = e lx ljv with some x G M. It 
implies that (rjLiVR^O-)* = (vl^ fjn)U (1)* , hence the verification of the injectivity of is complete. 

Finally, we are showing that is an immersion. For, let x = (X,9, (r]L,r] R )U(l)*) G Ai R be an 
arbitrary point and take some tangent vector 

v = SX © 69 © (vl,Vr)(X l © Xfl)u(l), G T X M R (4.47) 

generated by the tangent vectors 5X G T\t = R n , 59 G T e R n = R n and X L ,X R e t Our first goal is 
to find the action of the derivative 

(dT R ) x : T X M R -> T t h {x) &o (4.48) 
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on vector v. As in ()4.2ip . let 1Z denote the square root of A. Also, let us introduce the shorthand 
notations 5C = C(6X) £ p, p = £(V(A, 6)) € O, together with the tangent vectors 

51Z = (d1Z) (m (5X © 56) 6 T n{Xfi) G, (4.49) 
5T = (dJ) {xfi) (5X © 59) € T T[xfi) C N - C", (4.50) 

and 



5p=(d(^oV)\ m (5X®59) 

=\gll- 1 ((5T)T* + T(5T)* - (511)11^ TT* - TT*TV X (5V,)) 1Z' 1 e T p O. 



(4.51) 



Utilizing the above objects, one can easily verify that 

(dT R ) x v = rj L (5K + X L K - KXr)^ 1 © rj R (5C - [CXr])^ 1 © r, L (5p - [p, Xl})^ 1 . (4.52) 

To find the kernel of the linear transformation (dT R ) x , let us suppose that v € ker((dT R ) x ). 
Note that both C and 5C are diagonal matrices, meanwhile the commutator [C, Xr] is off-diagonal. 
Therefore, the vanishing of the g-component of (|4.52|) entails the relations 5C = and [£, Xr] = 0. 
Whence 5X = 0, and by the regularity of C the matrix Xr must be diagonal, i.e. Xr = diag(i%, i%) £ m 
with some x € Combining this fact with the vanishing of the G-component of (|4.52p . we obtain 

X L = TZXrU- 1 - {5U)n- 1 . (4.53) 

For simplicity, let us introduce the notation 

5 A = (dA\ m (8X © 56) e T A(Xfi) G. (4.54) 

Since A = 1Z 2 , the relationship between their derivatives leads to the formula 

5A = {511)11 + 11(511). (4.55) 

Keeping in mind (|4.5ip . ()4.53p and (|4.55p . the vanishing of the O-component of the tangent vector 
(|4.52|) immediately yields 

(5T)T* + T{5T)* - TT*X R + AXrA^TT* - (5A)A' 1 TT* = 0. (4.56) 

However, since 5X = 0, the derivative of the vector-valued function T fj4.9|) comes easily. Namely, for 
the components of 5T we simply have 

(5T) a = (56) a T a and (5F) n+a = -(56) a F n+a (o € N„). (4.57) 

Upon introducing the N x N diagonal matrix 

D = diag ((59)x, . . . , (56) n , -(56)^. . . , -(66) n ) e a, (4.58) 

the above relations (|4.57p can be cast into the simpler matrix form 5J- = DT . Moreover, by inspecting 
the derivative of A (|4.10p . we obtain the concise expression 

(5A)A~ X =D + ADA' 1 . (4.59) 
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Plugging these formulae into (|4.56[> . we get 

[X R - D, A' 1 TJ 7 *] = 0. (4.60) 

Remembering (|4.27p . we can write AC J- = —J 7 , whence A^ 1 !- = —CT. Therefore the commutation 
relation (|4,60p can be rewritten as 

(X R - D) k>k {CF) k Fi = (CJF) k Ti(XR - D\i (Vfc, I G Nat). (4.61) 
Since the components of J- are non-zero, it follows immediately that 

(X R - D) kyk = (Xr - £>),,, (Vfc, i G Njv). (4.62) 
In particular, for each a G N n we have 

iXa - (W) a = (Xr - D) a , a = (X fi - D) n+a , n+a = i X a + (*0) a , (4.63) 

whence (<5#) a = 0, i.e. 56 = 0. Plugging this back into (|4.62p . we get Xr = ixljv with some \ G M. 
Utilizing (|4.53p . we obtain that Xl = ^X^-N, too, therefore X^(BXr G u(1)*. That is, each component 
of the vector v ([4.47P is zero, meaning that the kernel of the derivative (dT R ) x is trivial. ■ 

4.3 Canonical coordinates on the reduced phase space 

Now we are in a position to perform the second, concluding step of the symplectic reduction. Without 
repeating our discussion in the paragraph preceding theorem 5, let us notice that the base manifold 
of the trivial fiber bundle 

7T R :M R^ V R 7 (A, 9, (rjL, ^ij)C^(l)*) -> (A, 6) (4.64) 

provides an appropriate model for the reduced phase space V ext //o(K x K). In order to obtain the 
reduced symplectic structure, it is tempting to imitate the proof of theorem 5. Since we had very 
explicit and simple formulae for the objects appearing in the definition of T s (|3.25|) . in that case the 
computation of the pull-back (Y s )*u; ext was almost trivial, hence by ()3.53p the reduced symplectic 
form was also immediate. However, the definition of (|4.38|) involves the square root of A, for 
which we have no explicit formula. Therefore a direct computation of the pull-back (T R )*uj cxt seems 
to be hopeless. To circumvent this apparent difficulty, we instead switch to the analysis of the reduced 
Poisson bracket {• , - } R induced by the reduced symplectic form u R via a formula analogous to ()2.34p . 

In what follows, we introduce two families of auxiliary smooth functions defined on the extended 
phase space J> ext (j2.28p . Namely, for each r G N we define 

^:P ext ^M, (y,Y,p) <p r (y,Y,p) = + tl ^ Y *Y) (4 g5) 

2r 

* r: ^ R , (y,y, P ) » * r[y ,Y, P ) = «T*zm + tK(ny*Md (466) 

where Z(p) = (ig)^ 1 p+ In — is a Hermitian N x N matrix. Take an arbitrary point x = (y, Y, p) G 
■p ext and an arbitrary tangent vector Sy © 8Y © 5p G T x V ext . Utilizing the bilinear form (|2.4p . for the 
derivative of ip r we can write 

if v is odd 

(difrUSy (B5Y® Sp) =5y(B6Y(B 5p[p r ] = { "l r _, „ rx . c . ' (4.67) 

it r is even. 
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A slightly longer computation also reveals that the derivative of ^f r has the form 

/ y t + (Y*Y Y r + (Y*Y \ 

(d*rM6y 5p) = l—L±—L y *Z{p)y - Cy*Z(p)y ^ ' C, y'Hyj 

| Y r - 1 -iy*Z(p)yYi - C{Y*y-^y*Z{p)y(Y*yC \ ^ 

+ / ypr + (Y*y)y* + c y (Y r + (Y*y)y*c § \ 

Recalling the definitions (|2.10p . (|2.27p . (|2.29)) and (|2.33j) . we can even determine the Hamiltonian 
vector fields X Vr G X(V ext ) and Xq, r E X('P ext ) generated by the functions ip r and \lv, respectively. 
Indeed, at the point x we have 



(U = iLu«.« ', (4-69) 



0, if r is odd, 

yy r_1 ©0®0, if r is even. 

It is also immediate that we can write 

(X 9r ) x = Ay (BAY® Ape T x V cxt (4.70) 

with components 
r-l 

A V = | E (^-V^Cp)^ - C(y*) r - 1 -^*Z(p)y(y*)^C) € T,G, (4.71) 
j=0 

AY = ± (C(Y*) r y*Z(p)yC + Cy*Z(p)yY r C - (Y*) r y*Z(p)y - y*Z(p)yY r ) G T Y g - g, (4.72) 

Ap = ^-[y(Y r + (yT)y* + Cy(Y r + (Y*Y)y*C,p] € T p O. (4.73) 
4i# 

To proceed further, we continue with some standard facts on the family of the if x if -invariant 
smooth functions 

C^ xK (V cxt ) = {He C°°(P ext ) | H o = for all (A*, e if x if}, (4.74) 

where <3? ext (|2.30|) stands for the natural action of the Lie group if x if on the extended phase space 
"P ext . As is well-known (see e.g. Theorem 4.3.5 in [E]), for each H € C^xicC^ 3 ) there is a unique 
H R G c°°(V R ) such that 

(tt r )*H r = (T R )*H. (4.75) 

The above function i7 R is called the reduced Hamiltonian associated with ff. In particular, the 
functions ip r ()4.65p and ()4.66p are obviously if x if -invariant. Making use of ()4.38p and ()4.64p . 
one can easily verify that the corresponding reduced Hamiltonians have the form 

R 1 0, if r is odd, 

bL a =i A a. if r is even, 



and 



= | 2 ELi A aM sinh(2fl a ), if r is odd, 
l 2 Ea=i K\ z a\ cosh(26» a ), if r is even. 
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From the perspective of symplectic reduction it is a crucial fact that the family of the K x i^-invariant 
smooth functions forms a Poisson subalgebra of C°°('P cxt ). Moreover, for any H\,H2 G C^ x ^(P ext ) 
we have 

{ir R )*{H R , H R } R = {T R )*{H 1 ,H 2 } cxt . (4.78) 
It is also immediate that the map 



Cf xf (r') 3H^H R G C°°(V R ^ 



(4.79) 



is a Poisson algebra homomorphism. Based on the above observations, in a sequence of propositions 
we will show that the global coordinates A a and 9 a (a G N n ) form a Darboux system on the reduced 
phase space V R (|1.3p . 

Proposition 11. For any a, b G N n we have {X a , \b} R = 0. 

Proof. Let r, s G N be arbitrary even numbers. Also, let x = (y, Y, p) G V be an arbitrary point; 
then by (j!T34"j) . (^67]) and (j4T69|) we can write 



{ip r , cp s } ext (x) = (x v ,) x [<p r ] = yY 8 - 1 e o e o[<^ r ] = o. 

Due to ()4.78p we obtain immediately that 

{ip R ,ip R } R = 0. 



(4.80) 



(4.81) 



However, the above bracket can also be analyzed by exploiting the bilinearity, antisymmetry, Jacobi 
identity and Leibniz rule, characterizing any Poisson bracket. Indeed, recalling (|4.76p . we can write 



VPr > <Ps I 



rs 



a,b=l a, 6=1 

Now let us introduce the Vandermonde-type n x n matrix QJ with entries 

^a,b = X 2 a b - 1 (o,6GN„). 
Since A G c, i.e. Ai > . . . > A n > 0, for the determinant of 5J we have 



(4.82) 



(4.83) 



det(QJ) 



Ai A^ ... A^ 
A 2 A2 ... A" 



2n-l 



2n-l 



A n A^ 



A 



2n-l 



n - a c 2 ) n Aa ^ o. 

l<c<d<n 



(4.84) 



0=1 



Also, we define the n x n matrix Vt with entries £l a b = {^a,Xf,} R (a, 6 G N ra ). Comparing equations 
(I4.8ip and (|4.82p . with the identifications r = 2c and s = 2d we can immediately write 



(4.85) 



a,6=l 



o,6=l 



o,6=l 



for all c,d E N n . These relations can be cast into the matrix form QJ*S123 = 0. However, from equation 
(|4.84p we see that QJ is invertible, therefore = 0, i.e. {A a , X b } R = for all a, b G N n . ■ 

Proposition 12. For any a, 6 G N n we /iai>e {^A;,}^ = \5 a ,b- 
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Proof. Take an odd number r and an even number s. Let x = (y,Y,p) G 'P ext be an arbitrary point 
with the additional assumption Y* = Y. Recalling (|2.34|) . (|4,68|) and (|4.69|) . one can easily verify that 

{* r , ^} ext (x) = (X v ,) x [*r] = VY 3 - 1 0[* r ] = 2* r+s _x(x). (4.86) 
Therefore, by (|4.78p . we obtain the relationship 

= 2*%..-v (4.87) 

From equations (|4.76p . (|4.77p . and from the previous proposition it is immediate that the left hand 
side of the above equation can be rewritten as 

a n n n 

{9?,<pf} R = -J2 Kka|sinh(20 a ),A^ = 4 £ X r a \z a \ cosh(20 a ) ^{20 a , X b } R X s b -\ (4.88) 

a, 6=1 a=l 6=1 

Comparing this expression with the right hand side of (|4.87|) . from (|4.77p we see that 

n / n \ 

^A^|z a |cosh(20 a ) ^{2^,A 6 } i? A^ 1 -Ar 1 =0. (4.89) 

o=l \6=1 / 

Upon introducing the n x n matrix H with entries 

S M = \z a \ cosh(20 a ) (fli 29 ^ h}**?' 1 ~ A^ 1 \ (a, d G N n ), (4.90) 

and recalling the Vandermonde-type matrix QJ introduced in ()4.83p . with the identifications r = 2c— 1 
and s = 2d the above equation (|4.89p takes the form 

n n 

= E A ' C " ls a,d = ^(2J*)c,aH a , d = (2TS) c , d , (4.91) 

a=l a=l 

for all c, d G N n . It simply means that = 0. However, due to (|4.84p the matrix QJ is invertible, 
hence 5 = also holds. Upon introducing the nxn matrix Q, with entries Jl a & = {29 a , X b } R (a, b G N n ), 
the equation z, a ,d = amounts to the requirement 

n n 

= J2{29 a , A^Af- 1 - A^- 1 = n a , b Wb,d - ^a,d = - 9J) 0)d , (4.92) 

6=1 6=1 

for all o,rf£ N n . In matrix notation we can simply write f2QJ = QJ. Again, since the Vandermonde-type 
matrix 23 is invertible, it entails Q = l n , i.e. {29 a , Xb} R = 8 a j for all a, 6 G N n . ■ 

Proposition 13. For any a, 6 G N n we /iai>e {0a>06} fi = 0. 

Proof. Though the proof of this proposition is computationally more demanding, the idea is the 
same as in the previous two propositions. Namely, take an arbitrary pair of odd numbers r, s £ N and 
compute the reduced Poisson bracket {^f R ,^ R } R in two different manners. 

First, let (A,0) G V R and define x = (A{X, 0)3, £(A), £(V(A, 9))) G P oxt . Recalling (IQ81 and 
(jITMD . by g^SD and (fX3H) we can write 

{^,*f} fl (A,0) = {Vr^^ix) = (Xv.) x [*r] = {d* r ) x (Xv.) x . (4.93) 
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However, recalling (|4.70j) . we have 

(X* s ) x = Ay © AY Ap, (4.94) 

where the components can be determined from (|4.7ip . f|4.72|> and ()4.73p . respectively. Now, by applying 
(|4.68|) on the components of the tangent vector (Xqt a ) x , we obtain 

n s—1 

AyffiOffiO[* r ] = ^ ^A^'Ar W ka(A)|| 26 (A)|( e 2 ^-(-iye- 2 ^)(e 2 ^ + (-iy e - 2 ^) (4.95) 

a,b=l j=0 

and the similar expression 

n r— 1 

© Ay © 0[* r ] = - S A a +iA r 1_i k«( A )lk6( A )K e2ea - (-l) j e- 2 ^)(e 2eb + (-l) J e - 2ei >), (4.96) 
a,6=l j=0 

together with the slightly more complicated derivative 

OffiO© A r* 1 - 4 V \ r \ s \z (X)\\zJX)\ ( ^a-^)^H^a + 29 b ) (A, + A») sinh(2fl a - 2g 6 ) \ 
0©0©Ap[* r j -4^A a A 6 K(A)||^(A)| ^ %2 + (Aa _ A&)2 Ag2 + (Aa + Xh)2 ) ■ 

(4.97) 

Plugging these formulae into ()4.93p . simple algebraic manipulations lead to the expression 

n 

¥? } R = - 2(r - a) £ A^-V«[ 2 sinh(40 a ) 

a=l 

2 y> / A^Ag|z a ||^| sinh(2^ a - 2^) _ X r a X s b \z a \ \z b \ sinh(2fl a + 26 b ) \ ( 498 ) 
9 afc ^ V(% 2 + (A a + X b f)(X a + A 6 ) (4 5 2 + ( Aa - A b ) 2 )(A a - X b ) J ' 

(^46) 

On the other hand, we can also work out {^ R ,^ R } R relying only on the defining properties 
(bilinearity, antisymmetry, etc.) of the Poisson bracket. Indeed, utilizing (|4.77p and the previous two 
propositions, one can easily verify that 

=4 £ (\Z\z a \ cosh(2fl a ) sinh(2fl b ) ^bN) _ A g smh(2fl a ) cosh(2fl & ) ^ 

a,6=l ^ a 6 ^ 

n 

+ 16 ^ A^A^|z a ||z fe |cosh(20 a )cosh(2^){^,^} i? . (4.99) 



a,b=l 



At this point we need the partial derivatives of \z a \ (a £ N n ). However, recalling ()4.ip . it is immediate 
that for any a,b € N n , a ^ b, we have 

1 0|z 6 | _ l<91n(|z b | 2 ) _ 4 ff 2 1 4 5 2 1 



l^foi 9A a 2 9A a 4 9 2 + (A a - A b ) 2 A a - A 6 4<? 2 + (A a + X b ) 2 X a + A 6 " 

Along the same line we obtain 

kl dA a ^ + A 2 A a V% 2 + (Aa-A fe ) 2 A a -Afe 4 ff 2 + (A a + A fc ) 2 A a + Aj • l ' J 
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Now plugging these derivatives into (|4.99p . we see at once that the resulting formula coincides with 
the right hand side of (|4.98|) . up to some terms involving the desired Poisson brackets {9 a , 9 b } R . More 
precisely, the comparison of (|4.99p and f|4.98|> yields 

n 

\ r a \ s b \z a \\z b \cosh{29 a )cosh(29 b ){9 a ,9 b } R = 0. (4.102) 

a,b=l 

In order to infer {9 a ,9 b } R from the above equation, let us introduce the n x n matrix Vl with entries 

n atb =\z a \\z b \cosh{29 a )cosh{29 b ){9 a ,9 b } R (a,b€N n ). (4.103) 

Recalling the Vandermode-type matrix QJ defined in ()4.83p , with the identifications r = 2c — 1 and 
s = 2d — 1 we can write 

n n 

o=Y, ^, c n a ,m, d = (w*)c,an a , b Wb,d = (2TW) c ,d, (4.104) 

a, 6=1 a, 6=1 

for all c, d € N n , i.e. we have the matrix equation 23*^203 — 0. By the determinant formula (|4.84p the 
matrix 23 is invertible, whence ft = 0. Therefore, due to (|4.103p . we end up with {9 a , 9 b } R = for all 
a,6eN n . ■ 
Putting together the contents of the previous three propositions we arrive at the most important 
technical result of the paper. 

Theorem 14. Up to an obvious rescaling, the globally defined functions X a , 9 a (a 6 N n ) provide a 
canonical coordinate system on the reduced symplectic manifold V R . More precisely, for the reduced 
symplectic form uj r £ fl 2 (V R ) we have ui R = 2 Yla=i d#a A dA a . 

Before turning to the consequences of the above theorem, we conclude this section with some remarks. 
First, recall that the existing proofs of the analogous theorem in the A n case are quite involved. 
They rely either on non-trivial analytic facts (see e.g. the proof of Theorem CI in [2]), or on equally 
non-trivial functional identities (see the proof of Theorem 2 in [5]). On the other hand, the applied 
techniques in proving propositions 11, 12 and 13 are of purely algebraic nature. We believe that 
our proofs can easily be adapted to the A n root system, providing an elementary and conceptually 
simpler understanding even in the A n case. However, we find it a much more appealing aspect of 
our approach that it can be naturally extended to understand the duality between the hyperbolic 
Sutherland and the rational RSvD models associated with the BC n root system. Since the BC n case 
is computationally even more demanding that the C n case, we wish to come back to the BC n problem 
in a later publication. 

5 Action-angle duality and its consequences 

Based on two different diagonalization procedures, in the previous two sections we worked out the 
symplectic reduction of the extended phase space "P cxt = G x g x O (|2.28p at the zero value of the 
momentum map J ext (|2.3ip . The derivation of the phase space of the Sutherland models relies on 
the KAK decomposition of the G-component of the level space £o = (J ext ) _1 ({0}), meanwhile the 
phase space of the RSvD model comes from the diagonalization of the g-component. Although the 
approaches are apparently different, still, the resulting symplectic manifolds (V s , u s ) and (V r ,oj r ) 
are two models of the same reduced symplectic manifold V ext //o(K x K). It is therefore immediate 
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that there is a natural symplectomorphism S: V s — > V R between V and V R . More precisely, S is a 
symplectomorphism making the diagram 



ext. 



V 



(5.1) 



V 



■v 



R 



commutative. Moreover, as we see in theorems 5 and 14, the standard coordinates q a , p a (a G N n ) 
and A a , 9 a (a G N n ) provide Darboux systems on V and V R , respectively. More precisely, they are 
canonical only up to some trivial rescaling. Nevertheless, by slightly abusing the terminology, in the 
rest of the section we shall simply refer to them as canonical coordinates. Now, for their pull-backs 
we introduce the shorthand notations 



A a — <S*X a , 



S*9 a , 



(s- 



Qa, Pa 



(<s- 



Pa, 



(5.2) 



for all a G N ra . Since S is a symplectomorphism, it is obvious that the functions X a ,0 a £ C°°(V S ) 
(a G N„) provide a new family of (global) canonical coordinates on V s , meanwhile the functions 
Qa, Pa £ C°°('P^) (a £ N n ) are new canonical coordinates on V R . At this point we are in a position to 
harvest some of the immediate consequences of the dual reduction picture ()5.ip . 

Starting with the Sutherland side of the dual reduction picture, take an arbitrary (q,p) G V s and 
consider the point 

x = T s (q,p, (1 N , l N )U{l) m ) = (e Q ,L(q,p),aE)) G £ - (5.3) 

Making use of the parametrization ()4.38p and introducing the shorthand notation (A, 0) = S(q,p) G 
V R , from (|5.ip it is immediate that 

x = (r /L ^(A,e)^ J j 1 ,7 ?i? £(A)^ 1 ,7 ?L e(V(A,^))7 ?L 1 ), (5.4) 

with some r]L,f]R G K. Comparing (|5.3p and (|5,4p it is clear that 

L(Q,P) = VR^WVr 1 , (5-5) 
therefore the spectrum of L(q,p) can be identified as 

a(L(q,p)) = {±X a (q,p)\aen n }. (5.6) 
Now take an arbitrary real- valued Ad-invariant smooth function F : g — >• M, i.e. we require that 

FiyYy- 1 ) = F(Y) (Vy G G, VY G g). (5.7) 
Recalling the notations introduced in (|5,2p . from the relationship (|5.5p we see that 

FoL = F(£(Ai,...,A n )), (5.8) 



i.e. the composite function F o L G C°°{V S ) depends only on the functions A a (a G N n ). Therefore 
the canonical coordinates A a , 6 a (a G N n ) provide an action- angle system for the Hamiltonian system 
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(V , oj , FoL). It is worth noting that the action variables A a are actually the positive eigenvalues of L 
(see (|5.6|) ). meanwhile the angle coordinates are the pull-backs of the impulses 9 a of the RSvD picture. 
It is a trivial consequence of the above observations that during the evolution of the dynamics governed 
by the Hamiltonian FoL, the Hermitian matrix L (|3,13p undergoes an isospectral deformation. 
Moreover, by (|5.6p . the positive eigenvalues of L provide n functionally independent smooth functions 
in involution, therefore L is indeed a Lax matrix for the system (V , ui s , FoL). Notice also that with 
the aid of the real-valued Ad-invariant function 

F 2 (Y) = \(Y,Y) = \tr{Y 2 ) (Y G q) (5.9) 

we can recover the hyperbolic C n Sutherland system with two independent coupling constants. Indeed, 
one can easily verify that 

F 2 oL = H% n , (5.10) 

therefore we have complete control over the Lax matrix and the action-angle variables of the model 
(jl.2p . We mention in passing that though the Lax matrix of the hyperbolic C n Sutherland model has 
been known for decades (see e.g. the survey [9|), to our knowledge the construction of action-angle 
coordinates for the non-^4 n -type Sutherland models has not been carried out in the literature before. 

Having finished our discussion on the Sutherland picture, we now develop a parallel theory in the 
dual RSvD picture, too. For, take an arbitrary (A, 9) G V R and let 

x = (A(A,0)3,/:(A),£(V(A,0))) G £ . (5.11) 

Recalling T s f|3.25[) and defining (q,p) = S~ 1 (X,9) G V s , from (|5.ip it is evident that 

x = (r] L e Q r] R ^ : r]RL(q : p)r] R ^,r] L ^(E)r]l 1 ) (5.12) 

with some i^LiVR € K. Comparing the G-components of (|5.11|) and (|5.12p . we obtain 

A(X,9) l 2 = r, L e^\ (5-13) 

which in turn immediately yields -4.(A, 9) = rjLe 2 ®^ 1 . Therefore the spectrum of the positive definite 
matrix A(\,9) has the form 

a(A(X, 9)) = { e ±2 ^ x ^ | a G N n }. (5.14) 

Now take an arbitrary real- valued K x i^-invariant smooth function /: G — > R. More precisely, we 
require the invariance of / under the K x i^-action (|2.15p . i.e. we impose 

f((k L ,k R ).y) = f(k L yk R 1 ) = f(y) (V(k L ,k R ) G K x K,Vy G G). (5.15) 

From (|5.13p it is obvious that 

foA* = /(diag(e ft ,...,e**,e _ *,...,e _ **)), (5.16) 

i.e. the composite function / o A* G C°°(V R ) depends only on the coordinates q a (a G N n ). It is 
now immediate that the canonical coordinates q a , p a (a G N n ) give rise to an action-angle system 
for the mechanical system (V R ,uj R ,f o A^). For purpose of interpretation we note that the action 
variables q a are provided by the positive eigenvalues of ln(^4)/2 (see (|5.14p ). meanwhile the angle 
variables come from the pull-backs of the canonical impulses p a appearing in the Sutherland picture. 
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Rephrasing the above observations, we see that the positive eigenvalues of ]n(A)/2 form a commuting 
family of n functionally independent smooth functions. Moreover, if the dynamics is generated by the 
Hamiltonian / o A? induced by some K x JT-invariant function /, then the matrix A (|4.2p undergoes 
an isospectral deformation. Therefore it is fully justified to call A a Lax matrix for the Hamiltonian 
system (V R ,u R ,f o A?). To make contact with the RSvD models, let us notice that the real-valued 
K x K-invariant function 

fi(y) = ~tr(yy*) (y G G) (5.17) 
generates the rational C n RSvD model f j 1 . 4 [) with two independent coupling parameters, i.e. we have 

fiaA3=H$ n . (5.18) 

To sum up, besides the Lax matrix A, from the dual reduction picture (|5.1|) naturally emerges an 
appropriate action-angle system of canonical coordinates for the most general rational C n RSvD model. 
Both the Lax matrix A and the construction of the action-angle variables for the C n RSvD system 
appear to be new results. 

Apart from a natural construction of action- angle coordinates, the symplectic reduction framework 
provides also nice solution algorithms for both the Sutherland and the RSvD models. Starting with 
the Sutherland picture, take an arbitrary Ad- invariant smooth function F: q —> R (|5.7p . Let 

pr g : V ext = Gxgx0-»g (5.19) 

denote the natural projection onto the factor g; then clearly pr*F = -Fopr G C^ xK {V ey± ). Notice that 
in the Sutherland picture the reduced Hamiltonian induced by the K x if-invariant smooth function 
pr*F is exactly FoL. Thus, as is known from the theory of the symplectic reductions (see e.g. Theorem 
4.3.5 in pi]), the Hamiltonian flows of pr*F flowing on the level space (■7 tnrt ) -1 ({0}) = £o = M s 
project onto the Hamiltonian flows of the reduced Hamiltonian system (V s , uj s , FoL). To understand 
the Hamiltonian flows of pr*F we need the g-valued gradient VF of F, which is defined by the formula 

(dF) Y (5Y) = (VF(Y),5Y) (Y G g,5F G TyQ = g). (5.20) 

Let us notice that the infinitesimal version of the Ad-invariance of F takes the form 

[Y,VF(Y)}=0 (VFGg), (5.21) 

whence from equations (|2.33p and (|2.29p we see that at each point (y, Y, p) £ P ext the Hamiltonian 
vector field X pi *p E X(7 ,cxt ) has the form 

(X pi * F ) (yXp) = (yVF(Y)) y 0y P G T {yXp) V cxt . (5.22) 

It is now immediate that the Hamiltonian flows of prg-F are complete and they have the form 

R3t^ (y e™ F{Yo \Y , Po ) G P ex \ (5.23) 

with some (yo;^0)Po) G p 6 ^. Therefore the flows generated by the reduced Hamiltonian FoL are 
also complete. Moreover, if 

R 3 t i-> (q(t),p(t)) G V S (5.24) 
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is an arbitrary flow of the reduced Hamiltonian F o L, and if we introduce the shorthand notation 
Lq = L(q(0),p(0)), then the flow of pr*F passing through the point (e™ \ Lq, £(E)) € £o projects 
onto the flow ()5.24p . More precisely, for each t 6 1 there is some (j]L(t),r]ii(t)) € K x K such that 

{e^e tWF ^\L^E)) = { nL {t)e^ (5.25) 

Let us observe that the g-component of the above equation has the form 

L = m (t)L(q(t),p(t)) m (t)-\ (5.26) 

which confirms the fact that during the time evolution of the dynamics the Lax matrix L (|3.13p does 
undergo an isospectral deformation. More interestingly, from the G-component of (|5.25p it is clear 
that 

e Q(0) e WF(L ) e tVF{L r e Q(0) = rj L ^) e 2 ^ VL (t)~\ (5.27) 

which immediately yields the spectral identification 

a(e 2Q V) = a (e 2Q ^e tyF ^e WF ^). (5.28) 

Our conclusion is that the diagonal matrix Q(t), and so the trajectory q(t) 6 c (t 6 R), can easily be 
recovered simply by diagonalizing the matrix flow 

t ^ e 2Q(0) e tVF(L ) e tVF(L r_ ( 5 . 29 ) 

In particular, by choosing the quadratic function F = F<i (|5.9p . we can reconstruct the flows of the 
hyperbolic C n Sutherland model (|1.2p by diagonalizing the exponential matrix flow 



3 2Q(0) JL 



(5.30) 



Here we used the obvious fact that Vi*2(Lo) = Lq/2 € p. 

Turning to the RSvD side of the dual reduction picture (|5.ip . take an arbitrary K x K- invariant 
smooth function f:G—>M. ()5.15p defined on the group G. If 

pr G : V cxt = GxgxO^G (5.31) 

denotes the natural projection onto G, then clearly pr G / = / o pr G € C-^ xK (V ext ), whence pr^/ 
survives the reduction. More precisely, in the Ruijsenaars picture the corresponding reduced Hamil- 
tonian turns out to be / o A? . Therefore the Hamiltonian flows of pr^/ staying on the level space 
(J cxt ) _1 ({0}) = £o — M R project onto the Hamiltonian flows of the reduced system (V r ,uj r , f oA^). 
To ease the calculations, let us introduce the g-valued gradient V/ of / by the requirement 

(df)y(Sy) = {^f(y),y~ 1 Sy) (y &G,5y€ T y G). (5.32) 

Utilizing the gradient, for the Hamiltonian vector field X pr * a f 6 X("P cxt ) at each point (y,Y,p) £ V ext 
we can write 

(X pihf ) {yX>p) =0 y (B (-V/(y))y © Op € T (yXp) V ext . (5.33) 
Clearly the flows generated by the Hamiltonian pr^,/ are complete, having the form 

R3t^ (y ,Y -tVf(y ),p ) € P cxt , (5.34) 
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with some (jjo, Yq, pq) £ -p ext . Therefore the flows of the reduced Hamiltonian foA 2 are also complete. 
Now take an arbitrary Hamiltonian flow 

R3t^ (\(t),6(t)) £V R (5.35) 

of the reduced system (V R ,uj R ,f o A 2 ). By realizing (|5.35[) as a projection of an appropriate flow 
(|5.34p on the levels space £q = M. R , \t is clear that for each tslwe have 

(4,A-tv/(4u(v )) 

= (r ?L (t)^(A(t),^t))l ??R (0- 1 ,^(t)/:(A(t)) % (t)- 1 , ??L (t)e(V(A(t)^(t)))r /L (0- 1 ) (5.36) 

with some r) L {t),r) R (t) € K, where A = A(X(0), 0(0)), £ = £( A ( )) and V = V(A(0), 0(0)). By 
inspecting the G-component of (|5.36p it is immediate that 

A = VL(t)A(X(t), 9{t))r, L (tY\ (5.37) 

which clearly shows that the evolution of the Lax matrix A (|4.2p is isospectral. Prom our perspective 
it is more interesting that the g-component of (|5.36p has the form 

Co ~ tVf(Al) = rjR^CiXit^rjRit)- 1 , (5.38) 
from where we obtain the spectral identification 

«(C(\(t))) = a(C - iV/(4))- (5-39) 

The outcome of our analysis is that the diagonal matrix £(A(t)), and so A(t) € c (t S M), can be 
determined by diagonalizing the linear matrix flow 

t^Co-tVf(Al). (5.40) 

In particular, with the choice f = fi (|5.17p the above diagonalizing procedure permits us to construct 
the trajectories of the most general rational C n RSvD model (|1.4|) . 

Based on the action-angle duality and the above simple solution algorithms one can build up the 
scattering theory for both the Sutherland and the RSvD models. In this respect the main thrust comes 
from the intimate relationship between the scattering properties of the Calogero-Moser-Sutherland- 
type many-particle systems and the scattering characteristics of certain integrable soliton equations. 
The connection between the ^4 ra -type particle systems and the soliton equations defined on the whole 
line is well understood (see e.g. [19], [20] and [21]), but the link between the non-A n -type particle 
systems and the soliton systems defined on the half-line is far less elaborated (see e.g. [22J). Motivated 
by this intriguing relationship, the study of the temporal asymptotics of the exponential matrix flow 
(|5.30|) in our recent paper [IB] led us to the conclusion that the hyperbolic C n Sutherland model is a 
pure soliton system with a factorized scattering map. In the same spirit, by exploring the temporal 
asymptotics of the matrix flow (|5.40p . one can work out the scattering map of the rational C n RSvD 
model as well. However, a decent analysis would require the study of the M0ller wave transformations, 
too, so to keep the paper in a reasonable length we defer the scattering theoretic analysis into a later 
publication. 
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To conclude this section, let us notice that apart from its interest in the theory of integrable 
many-particle systems and soliton equations, our results on the action-angle duality of the C n -type 
models may have applications in the theory of random matrices as well. Built on the Lax matrices 
of the Calogero-Moser-Sutherland-type many-particle systems associated with the A n root system, 
the authors of the recent papers [23] and [24] have constructed certain integrable random matrix 
ensembles characterized by spectral statistical properties 'intermediate' between the Poisson and the 
GUE statistics. It turns out that in their investigation the main technical tool is the action-angle 
duality between the Calogero-Moser-Sutherland and the Ruijsenaars-Schneider integrable models of 
type A n . Utilizing the Lax matrices L (|3,13p . A (|4.2|> . and the dual reduction picture (|5.ip . it appears 
to be a challeng ing exercise to extend their considerations onto root systems other that A n . 



A Cauchy matrices and functional identities 

To understand the symplectic aspects of the rational C n RSvD model, especially the properties of its 
Lax matrix (|4.2[) , we need some functional identities involving certain rational functions closely related 
to the theory of Cauchy matrices. The following identity, that we have learned from the appendix of 
[25], is of great importance in the study of both the Cauchy matrices and the RSvD models. 

Lemma Al. Let M, N € N be arbitrary positive integers satisfying M < N , and take two families of 
arbitrary complex numbers a±, . . . , cxm and . . . , (3n with j3 k ^ fii for k ^ I. Then for any complex 
number z G C \ {/3±, . . . , /3jv} we have 

nM I \ N n M in 

k=1 {z-a k ) _ ^ 1 U k =iW- a k 



>M,N 



uli(z - a) " JU,JV ' U * - Pi nil & - Pi) 



Notice that the above identity can be verified by simple residue calculus. 

Now, with any complex A^-tuple x = (xi, . . . , acjv) £ C N we associate the standard N x N Cauchy 
matrix C(x) with entries 

C k ,i(x) = — — (Me (A.2) 

1 + x k - xi 

The characteristic properties of the Cauchy matrices can be found in any advanced textbook on linear 
algebra (see e.g. [26J). Recall that the determinant of C(x) is given by 

det(C(x))= J] i-(J- a /)-2 - ( A - 3 ) 

l<fc<K7V 1 [Xk Xl) 

In order to get a concise formula for the inverse of the Cauchy matrix C(x), for each j £ Nat we define 
the rational function 

N 1 

= n -i^r ■ ( A - 4 ) 
k=i j k 

A direct application of (jA.lj) immediately leads to the identity 

.7=1 J 
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Now, let us introduce the N x N diagonal matrix 

W{x) = &\&g(w\(x), . . . , wn{x)). (A. 6) 

Utilizing (jA.ip . one can easily verify e.g. by induction on N that 

tr(W(a?)) = N. (A.7) 

Most importantly, with the aid of the diagonal matrix ()A.6P filled in with the rational functions (|A.4|) . 
the inverse of the Cauchy matrix C(x) can be computed by the formula 

C(x)- 1 = W{-x)C{-x)W{x). (A.8) 

By specializing the above identities, we can produce some particularly useful relations suitable for 
analyzing the Lax matrix of the C n RSvD model. 

Proposition A2. Let n £ N, N = 2n, and take an arbitrary complex vector x = (xi, . . . , xjv) G 
satisfying 

x n+a = -x a (Va G N n ). (A.9) 



Under this assumption the rational functions (A. 4) can be rewritten as 



w a (x) = w n+a (-x) = (l + TT ( 1 + ) i l + 7 ) ( a G N ™)> ( A - 10 ) 

V 2x aJ V; V Xa-X d J V X a +X d J 



d=l 
(<¥a) 



and for each k G Njv they enjoy the functional identities 

i \ N 



y^A x L = and y ^> ±_ (AU) 

Z ^l + 2x j ^ (1 + Xj -x k )(l + 2 Xj ) l-2x k 



Moreover, the diagonal matrix W(x) HA. 6\) and the Cauchy matrix C{x) HA.%\) generated by the complex 
vector x fA. 9\) satisfy the relationships 

CW{x)C = W(-x), CC(x)C = C(-x), {C{x)CW(x)f = 1 N . (A.12) 



Proof. We prove only the identities displayed in (jA.lip . With the identifications M = N , z = —1/2, 
a\ = xi — 1, fti = xi (I G Njv)j the application of the identity (jA.lj) immediately yields 

f '"jW _ !f 1 TT Jj-^-l) ,!^ Tr (l C -;)(ln+c-^ _ . 

Z. 1 + 2Xj ^^(-D-^.ll 2^ il(, c+ i) ( x n+c + i)J a (A - i3j 

Notice that the assumption (|A.9P on x is crucial in the last step of the algebraic manipulations. Finally, 
take an arbitrary number k G Njv; then from (|A. 13|) and (|A.5|) it follows that 

y w i( x ) _ 1 l 2 y ^iM w i( x ) 1 _ 1 ( AU) 

+ -x k ){l + 2x j ) l-2x k \^l + 2xj f^l + Xj-XkJ ~ l-2x k 1 ' 

The rest of the proposition is trivial. ■ 
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